A control volume finite element method for three-dimensional three-phase flows

A novel control volume finite element method with adaptive anisotropic unstructured meshes is presented for three-dimensional three-phase flows with interfacial tension. The numerical framework consists of a mixed control volume and finite element formulation with a new P₁DG-P₂ elements (linear discontinuous velocity between elements and quadratic continuous pressure between elements). A “volume of fluid” type method is used for the interface capturing, which is based on compressive control volume advection and second-order finite element methods. A force-balanced continuum surface force model is employed for the interfacial tension on unstructured meshes. The interfacial tension coefficient decomposition method is also used to deal with interfacial tension pairings between different phases. Numerical examples of benchmark tests and the dynamics of three-dimensional three-phase rising bubble, and droplet impact are presented. The results are compared with the analytical solutions and previously published experimental data, demonstrating the capability of the present method.     

On the inflation and deflation dynamics of liquid-filled,hyperelastic balloons

We derive a reduced-order model describing the inflation and deflation dynamics of a liquid-filled hyperelastic balloon, focusing on inviscid laminar flow and the extensional motion of the balloon. We initially study the flow and pressure fields for dictated motion of the solid, which throughout deflation are obtained by solving the potential problem. However, during inflation, flow separation creates a jet within the balloon, requiring a different approach. The analyses of both flow regimes lead to a simple piecewise model, describing the fluidic pressure during inflation and deflation, which is verified by finite element computations. We then use a variational approach to derive the equation describing the interaction between the extensional mode of the balloon and the entrapped fluid, yielding a nonlinear hybrid oscillator equation. Analytical and graphical investigations of the suggested model are presented, shedding light on its static and dynamic behaviour under different operating conditions. Our simplified model and its underlying assumptions are verified utilizing a fully coupled finite element scheme, showing excellent agreement.     

基于平面度的燃料电池电化学性能模拟

 固体氧化物燃料电池的翘曲会影响电极-盖板界面的接触情况，从而影响电化学性能，对相关制造工艺提出了很大的挑战．为了分析燃料电池平面度对放电过程的影响，揭示其潜在的风险，我们建立了两个基于有限元法的仿真模型，对考虑平面度缺陷的燃料电池封装和放电进行分析．在对固体氧化物燃料电池进行平面度测量的基础上，首先建立了具有真实燃料电池翘曲特性的几何模型，分析封装过程中接触压力的分布情况．然后将接触电阻的仿真结果导入到三维多物理场耦合模型中，模拟具有平面度缺陷的燃料电池电化学性能．计算结果展示了燃料电池两侧封装过程中接触压力的分布情况．通过对比有接触电阻和无接触电阻的燃料电池电流密度，分析了电池与盖板的接触对放电过程的影响．结果表明，燃料电池的凹陷面较难达到满意的接触状态，需要比凸起面更大的封装压力．燃料电池表面接触电阻的变化将引起电流传导路径的变化，产生局部高电流或低电流．这项工作强调了在燃料电池中保持均匀分布的接触电阻的重要性，为今后的优化工作奠定了基础．     

Stochastic sensitivity analysis of coupled acoustic-structural systems under non-stationary random excitations

This paper presents a new sensitivity analysis method for coupled acoustic–structural systems subjected to non-stationary random excitations. The integral of the response power spectrum density (PSD) of the coupled system is taken as the objective function. The thickness of each structural element is used as a design variable. A time-domain algorithm integrating the pseudo excitation method (PEM), direct differentiation method (DDM) and high precision direct (HPD) integration method is proposed for the sensitivity analysis of the objective function with respect to design variables. Firstly, the PEM is adopted to transform the sensitivity analysis under non-stationary random excitations into the sensitivity analysis under pseudo transient excitations. Then, the sensitivity analysis equation of the coupled system under pseudo transient excitations is derived based on the DDM. Moreover, the HPD integration method is used to efficiently solve the sensitivity analysis equation under pseudo transient excitations in a reduced-order modal space. Numerical examples are presented to demonstrate the validity of the proposed method.     

水流损失和热补偿共同作用对增强型地热系统(EGS)产能影响的研究

基于青海共和盆地-3705m地热田实测数据,结合流固耦合传热理论并运用Comsol软件,建立了离散型裂隙岩体流体传热模型。考虑水流损失和热补偿共同作用,模拟得到了开采过程中上、下岩层(盖层和垫层)为绝热不渗透、传热不渗透、渗透传热时,储层(上、下岩层和压裂层)温度场的变化特征,分析了产出流量、水流损失、产出温度、产热速率的变化规律。研究结果表明:采热过程中产出流量始终小于注入流量;产出流量增幅速率先增大后减小,最后趋于稳定,前3a产出流量增幅超过总增幅量的3/4;忽略水流损失,将高估产热速率,采热初期甚至达到考虑水流损失时产热速率的3倍以上;考虑水流损失,产热速率呈先快速上升再趋于稳定后逐渐下降的趋势,最优开采时间为3a^11a;研究上、下岩层对产出温度的影响,仅考虑传热,采热寿命延长5.43%,同时考虑渗流传热时,采热寿命延长2.71%;采热前9a,水流损失占主导作用,即流入上、下岩层水流损失对产热速率的影响高于热补偿效应,开采10a后,热补偿效应占主导作用;同时考虑水流损失和热补偿效应得到的产热速率变化规律与实际工程更为符合,建议选择低渗透能力的上、下岩层延长增强型地热系统(EGS)运行时间。     

Existence of multiple and sign-changing solutions for a second-order nonlinear functional difference equation with periodic coefficients

In the present paper, we apply the method of invariant sets of descending flow to establish a series of criteria to ensure that a second-order nonlinear functional difference equation with periodic boundary conditions possesses at least one trivial solution and three nontrivial solutions. These nontrivial solutions consist of sign-changing solutions, positive solutions and negative solutions. Moreover, as an application of our theoretical results, an example is elaborated. Our results generalize and improve some existing ones.     

Information percolation and cutoff for the random‐cluster model

We consider the random‐cluster model (RCM) on with parameters p∈(0,1) and q ≥ 1. This is a generalization of the standard bond percolation (with edges open independently with probability p) which is biased by a factor q raised to the number of connected components. We study the well‐known Fortuin‐Kasteleyn (FK)‐dynamics on this model where the update at an edge depends on the global geometry of the system unlike the Glauber heat‐bath dynamics for spin systems, and prove that for all small enough p (depending on the dimension) and any q>1, the FK‐dynamics exhibits the cutoff phenomenon at with a window size , where λ_∞ is the large n limit of the spectral gap of the process. Our proof extends the information percolation framework of Lubetzky and Sly to the RCM and also relies on the arguments of Blanca and Sinclair who proved a sharp mixing time bound for the planar version. A key aspect of our proof is the analysis of the effect of a sequence of dependent (across time) Bernoulli percolations extracted from the graphical construction of the dynamics, on how information propagates.     

梯级溃坝洪水洪峰增强机制

我国在多条河流上修建了大量梯级水库, 梯级坝溃决诱发洪水大大超过单坝溃决洪水洪峰, 因此亟需加深对梯级坝溃决洪水洪峰增强机制的认识. 本文建立了梯级坝溃决洪水演进过程的一维浅水动力学模型, 发展了一套能捕捉激波、干湿边界和保平衡结构的数值求解方法, 通过大量算例, 系统研究了梯级坝溃决洪水演进过程的质量转化和能量转化机制. 研究结果表明, 梯级溃决中, 上游溃决诱发的洪水大大增大下游水库的质量和动量, 形成一个带动量的水塔, 同时在尾部残留一个动量较大的射流, 不断补充下游坝体溃决后水塔的质量和动量, 持续维持洪峰高度. 根据该射流-水塔机制, 建立了梯级坝溃决洪水演进过程对应的射流-水塔单坝溃决洪水过程等效模型, 该等效模型基本反映了梯级坝溃决诱发洪水的洪峰过程, 并成功预测了多个坝间距为百公里量级的梯级坝溃决洪水洪峰高程和流量, 可望为流域防洪和梯级坝设计提供理论依据.